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Abstract 

An enhanced integrated aerodynamic load/dynamic 
optimization procedure, developed for minimizing 
vibratory root shear forces and moments of a helicopter 
rotor blade is described. The optimization problem is for- 
mulated with 4/rev vertical and 3/rev inplane shears at the 
blade root as objective functions. Constraints are imposed 
on 3/rev radial shear, 3/rev flapping and torsional mo- 
ments, 4/rev lagging moment, blade natural frequencies, 
weight, autorotational inertia, centrifugal stress and rotor 
thrust. The ’Global Criteria Approach’ is used for for- 
mulating the multi-objective optimization. Design vari- 
ables include spanwise distributions of blade bending 
stiffnesses, torsional stiffness, nonstructural mass, chord, 
radius of gyration and blade taper ratio. The program 
CAMRAD is coupled with an optimizer, which consists 
of the program CONMIN and an approximate analysis. 
The optimization procedure is applied to an advanced 
rotor as a reference design. Optimum blade designs, ob- 
tained with and without a constraint on the rotor thrust, 
are presented and are compared to the reference blade. 
Substantial reductions are obtained in the vibratory root 
forces and moments. As a byproduct, improvements are 
also found in some performance parameters which were 
not considered during the formulation of the optimization 
problem. The effect of thrust constraint on the values of 
the vibratory forces and moments is demonstrated by 
varying the magnitude of the prescribed thrust. A proper 
choice of the ’move limit’ paremeter, used in the ap- 
proximate analysis is shown to have significant effect on 
the optimum results. 
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autorotational inertia, lb-fr 

chord, ft 

root chord, ft 

tip chord, ft 

power coefficient 

thrust coefficient 

propulsive force coefficient 

bending stiffnesses, lb-fr 

bending stiffnesses at blade root, lb- ft “ 

natural frequencies of first four coupled 

elastic modes, per rev 

3/rev radial shear, lb 

3/rev inplane shear, lb 

4/rev vertical shear, lb 

objective functions 

approximate objective function 

inplane component of blade airload, lb/ft 

normal component of blade airload, Ib/ft 

f h constraint function 

approximate constraint function 

torsional stiffness, lb-ft 1 2 

torsional stiffness at blade root, lb-fr 

principle radius of gyration at blade root, ft 

maximum lift to drag ratio of an airfoil 

3/rev torsional moment, lb-ft 

3/rev flapping moment, lb-ft 

4/rev lagging moment, lb-ft 

move limit 

integer, 0 < n < 1 

number of blade nodes 
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NCON 

number of constraints 

NDV 

number of design variables 

NSEG 

number of blade segments 

R 

blade radius, ft 

Si 

centrifugal stress in i h segment, lb/ft 2 

Smax 

maximum allowable stress, lb/ft 2 

T 

thrust, lb 

Tref 

thrust of reference rotor, lb 

Wj 

nonstructural weight per unit length at j h 
node, Ib/ft 

W 

total blade weight, lb 

x, y, z 

reference axes 

a 

angle of attack, deg 

X 

taper ratio 

< Pi 

i th design variable 

4 

advance ratio 

a 

thrust- weighted solidity 

'P 

blade azimuth angle, deg 

a 

rotor angular velocity, rad/sec 


Introduction 

Vibration has been a major source of problems in 
helicopters and its alleviation plays an important role in 
the rotor blade design process. Due to more stringent re- 
quirements in the permissible vibration levels and require- 
ments for increased reliability* the helicopter industry is 
exploring methods to reduce vibration. In the past, con- 
ventional design methods mainly used the designer’s ex- 
perience and trial and error methods. Using these 
methods, rotor vibratory responses are reduced through 
post design addition of tuning masses, which incorporate 
weight penalty. Today with efficient optimization 
schemes available and improved helicopter analysis, at- 
tempts are being made to apply optimization techniques, 
at early stages of design, to the vibration problem. 

Due to the importance of the problem, there has been 
a considerable amount of research aimed at reducing 
vibration using optimization, primarily at the blade level, 
as shown in Refs. 1-11. In most of this work, the vibra- 
tion is controlled by reducing the oscillatory hub shears 
and moments. In one such effort, the mode shape is 
changed to make it orthogonal to the aerodynamic forcing 
function and thereby reduce vibration 1 . An early attempt 
at optimum blade design for reduced vibration is due to 
Peters et al. 2 where the objective was to place the natural 
frequencies away from the critical frequencies and there- 
by lower the hub shear. Refs. 3 and 4 present more 
recent developments in optimum blade designs for mini- 
mum weight with constraints on multiple natural frequen- 
cies, blade stress and autorotational inertia for articulated 
rotor blades with both rectangular and tapered planforms. 
However, the blade was assumed to be in vacuum and 


therefore the optimization was performed without air- 
loads. The results of the above research indicate that by 
appropriately selecting the blade inertial, structural and 
sometimes aerodynamic characteristics, it is possible to 
minimize blade vibration. 

In most of the previous work dealing with optimum 
blade designs for reduced vibration, the aerodynamic 
loads on the blade were prescribed and the effects of the 
design changes, during optimization, on the blade airloads 
were not included. In Ref. 10, for the first time, the in- 
tegration of aerodynamic loads and dynamics was ad- 
dressed by coupling a comprehensive helicopter analysis 
code, CAMRAD 11 , the nonlinear optimization algorithm, 
CONMIN 12 , and an approximate analysis technique. A 
combination of the blade root 4/rev vertical shear and the 
blade weight was minimized with constraints on coupled 
elastic flap/lag natural frequencies, blade autorotational 
inertia and centrifugal stress. The use of the program 
CAMRAD permitted the design of the blade with calcu- 
lated airloads and its presence in the closed-loop op- 
timization procedure allowed the inclusion of the effects 
due to changes in these airloads with changes in design 
variables. The paper demonstrated significant reductions 
in the 4/rev vertical shear and blade weight, which were 
objective functions. As a byproduct, it was shown that 
optimization also reduced the total power loading re- 
quired. 

The research reported in Ref. 15 extends the work of 
Ref. 10 by including other sources of blade vibration that 
are transmitted to the hub. The purpose of this paper is to 
include details which were not included in Ref. 15 due to 
page limitations. In Ref. 15 a more comprehensive for- 
mulation was presented by including the 3/rev inplane 
and radial shears, the 3/rev flapping and torsional mo- 
ments and the 4/rev lagging moment in the optimization 
formulation in the form of objective functions and/or con- 
straints. In Ref. 10, the optimum rotor was required to 
maintain the same Cr/cr and Cx/a as the reference rotor 
through a trim procedure, Ct being the rotor thrust coeffi- 
cient, Cx the propulsive force coefficient and a the thrust 
weighted solidity of the blade. This allowed the thrust of 
the optimum rotor to be lower than that of the reference 
rotor due to the smaller solidity of the optimum blade. 
To avoid this loss in the rotor thrust, in Ref. 15, an addi- 
tional constraint is imposed on the total thrust. 

Problem Description 

Optimization techniques are applied to minimize the 
vibratory blade loads of a four-bladed helicopter in for- 
ward flight. This is done by including the major sources 
of vibration (forces and moments) in the rotaung frame 
that are transmitted to the nonrotating system as 4/rev 
harmonics of the longitudinal, lateral and vertical forces 
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and the pitch, roll and yaw moments. For instance, in a 
four-bladed rotor, the 4/rev vertical shear at the blade root 
is transmitted as a 4/rev vertical force at rotor hub and the 
3 and the 5/rev components of the inpiane and the radial 
shears are transmitted as 4/rev forces in both lateral and 
longitudinal directions at the hub. Similarly, the 3 and 
the 5/rev components of the flapping and the torsional 
moments at the blade root affect the rotor hub in the form 
of oscillatory 4/rev pitch and roll moments. Finally the 
4/rev lagging moment at the root causes oscillatory 4/rev 
yawing moment at the hub. Therefore, in this study, all 
of the 3 and the 4/rev components of these critical shear 
forces and moments are included in the optimization for- 
mulation in the objective functions and/or constraints. 
Although the 5/rev components of the inplane shear, the 
flapping and the torsional moments are important, their 
values are much smaller and therefore are not included. 
The excessive addition of weight, associated with most 
vibration reduction processes, is avoided by incorporating 
an upper limit constraint on the blade weight. 

Optimization Formulation 

The objective functions are the 4/rev vertical shear, 
fz, and the 3/rev inplane shear, f x , at the blade root. The 
following constraints are imposed: 1) upper and lower 
bounds on first four elastic coupled blade natural frequen- 
cies, f3, f4, fs and f6 , 2) lower bound on blade autorota- 
tional inertia, AI, 3) upper bound on total blade weight, 
W, 4) upper bound on centrifugal stress on each blade 
segment, Si, i = 1,2,.., NSEG (NSEG is the number of 
blade segments), 5) upper bound on 3/rev radial shear, f r , 
6) upper bound on 3/rev flapping moment, m x , 7) upper 
bound on 3/rev torsional moment, m c , 8) upper bound on 
4/rev lagging moment, m z , 9) lower bound on total thrust, 
T and side constraints on design variables to avoid im- 
practical solutions. The bounds are selected from a 
baseline design which will hereafter be called the 
'reference* blade. It is required that the optimum blade 
maintains a certain level of the reference blade thrust. 
This leads to a lower bound constraint on the total thrust, 
T, which is expressed as 
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subject to the complete set of constraints 

gj ( 9 ) — 0 j=\X...,NCON 

and the side constraints. The design variable vector <pi * 
is obtained by minimizing the single objective function 
fz (cp) subject to the set of constraints g (rp) and the design 

variable vector <f >2 * is obtained by minimizing f x ((p) sub- 
ject to the same set of constraints. In the present case, 
due to the nonlinear nature of the above function the ob- 
jective function in Eqn. 2 is modified as follows : 

F((p) = ^F(q>) (3) 

The design variables are the blade stiffnesses at the root, 
EIxx n EIzz r > GJr, the taper ratio, X, the root chord, c r , the 
radius of gyration at the blade root, k r and the nonstruc- 
tural weights wj, j= 1,2,..., N, where N is the total number 
of blade nodes. 

Blade Model 

In this section the modeling assumptions that are used 
to relate the design variables to the blade structural 
properties like weight, autorotational inertia and 
centrifugal stress are described. A linear taper is assumed 
along the blade planform, as shown in Fig. 1 (to reduce 
the number of design variables) and the blade taper ratio 
is expressed as follows : 

(4 > 

Ct 

where c r is the root chord and Ct is the tip chord. It is 
assumed that the blade stiffnesses arise solely from the 
blade structural component and the contribution of the 
nonstructural masses, the skin and honeycomb are as- 
sumed to be negligible. The following radial distribu- 
tions are used for the stiffnesses: 
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fz ( 9 ) -fz (<pi) 

fz (cpi) 


T>nTref (o<n< 1) (1) 

where T re f is the reference blade thrust. A method called 
the ’Global Criteria approach* 13 is used to formulate the 
multiple objective optimization problem. Using this 
method, the optimum solution is obtained by minimizing 
a ’global criterion’ defined by the sum of the squares of 
the relative deviations of the individual objective func- 
tions from their respective individual optimum values. 
The optimization problem reduces to minimizing the 
global objective function, F (cp), where 
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where R is the blade radius and c(y) is the chord distribu- 
tion. For a linear taper 
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Similar distributions are used for EIzz and GJ. The link 
between the intermediate design variables and the blade 
cross sectional properties like area, weight , autorotational 
inertia and centrifugal stress are described in detail in 
Ref. 10. 

Analysis 

The program CAMRAD is used for both blade 
dynamic and aerodynamic analyses. The program CON- 
MIN, along with an approximate analysis technique, is 
used for optimization. 

Dynamic Analysis 

In CAMRAD the blade response is computed using 
rotating, free-vibration modes, equivalent to a Galerkin 
analysis. Ten bending modes, out of which seven are 
flapping (one rigid and six elastic) and three are lead-lag 
(one rigid and two elastic) and one rigid body torsion 
mode, are calculated. Blade resonances up to eight per 
revolution are included and therefore eight harmonics of 
the rotor revolution are retained in the air loads calcula- 
tion. The blade loads are generated using the isolated 
rotor model of CAMRAD. 

Aerodynamic Analysis 

CAMRAD is also used to predict the aerodynamic 
loading on the blade. Some assumptions made, while 
running the code, include: uniform inflow, yawed flow on 
the rotor, unsteady aerodynamics model and no dynamic 
stall. A trim analysis is performed inside CAMRAD at 
each step of design optimization. Since the reference 
blade is a wind tunnel blade model, one of the wind tun- 
nel options was selected for trimming the blade. The 
blade was trimmed to maintain the same Ct/cf, CJa and 
flapping angles as the reference blade using the collective 
and the cyclic (lateral and longitudinal) pitch and the 
shaft angle. 

Optimization Implementation 

The basic algorithm used is the method of Feasible 
Directions as implemented in the optimization program 
CONMIN. Since the optimization process requires many 
evaluations of the objective function and constraints 
before an optimum design is obtained, the process can be 
very expensive if full analyses are made for each function 
evaluation. The objective function and constraints are 
therefore approximated by linear Taylor series expansions 
based on the design variable values from CONMIN and 
the sensitivity information from the full analysis. Since 
the objective function and the constraints are all 
linearized, the optimization problem essentially reduces to 
a sequential linear programming (SLP) procedure. 
Specifically, if the objective function F, the constraint g, 


and their respective derivatives are calculated for the 
design variable cp* using an exact analysis, their values 
for an increment in the design variable A (p* are as fol- 
lows: 

NDV (7) 

^Acp* 


NDV 

k = 1 


where the quantities denoted ( A ) represent approximate 
values and NDV denotes the number of design variables. 
The assumption of linearity is valid over small increments 
in the design variable values and does not introduce large 
errors. A ‘move limit’ (ML), defined as the maximum 
fractional change of each design variable value, has been 
imposed as upper and lower bounds on each design vari- 
able cp*. In this paper, equal values are used for both 
upper and lower bounds. 

Test Problem 

The reference blade chosen for this work is a 
modification of a wind tunnel model of an advanced ar- 
ticulated rotor blade of the Growth Utility Rotor 14 . The 
modifications were made to the model to make it more 
suitable for optimization studies. For example, a single 
airfoil is used throughout the blade span instead of the 
three different airfoils present in the advanced blade. 
Also, a rectangular planform is assumed for the reference 
blade rather than the tapered. The reference blade has a 
radius, R = 4.685 ft, rotational speed, £2 = 639.5 rad/sec, 
flap hinge offset value (normalized with respect to the 
radius) of 0.0534 and maximum linear twist of -16 
degrees. Further details of the blade model can be found 
in Ref. 10. Optimum designs are obtained for a forward 
flight case with an advance ratio, ja = 0.3. 

Results 

In this section, results obtained by using the optimiza- 
tion procedure are presented. Optimum designs are ob- 
tained in 10-15 cycles. Optimum blade designs are ob- 
tained with a thrust constraint of T> T re f , he., the op- 
timized blade must maintain at least the reference blade 
thrust (n = 1.0 in Eqn. 1). Note that when T = T rc f, the 
solidities of the optimum and the reference blades will be 
the same since Ct/ct and Cx/cr are maintained constant 
through the trim procedure. Results are also obtained 
without the thrust constraint This represents a case 
where the thrust can be reduced if a is reduced, as long 


4 



as the quantities Ct/g and Cx/cr are maintained constant 
(in this case through the trim procedure), i.e., the thrust 
and the propulsive force per unit area are held constant 
during the optimization. The remaining constraints are 
the same for both cases. In addition results are also 
presented with varying magnitude of the required thrust, 
i.e., for different values of n in Eqn. 1. For ease of refer- 
ence, the following notation will be used 

Case 1: optimum design with thrust constraint, 
T > Tref 

Case 2: optimum design without thrust constraint. 

Clearly Case 1 is of much more practical interest. How- 
ever, Case 2 is included because it brings out some inter- 
esting phenomenological effects which illustrate impor- 
tant aspects of aerodynamic/dynamic interaction in 
rotorcraf. These integrated aspects are important to 
predict and understand in future efforts in integrated 
rotorcraft analysis and optimization. 

Optimization Results 

Table 1 presents a summary of the reference and the 
optimum designs in terms of the objective functions and 
constraints. By definition, a constraint is active when its 
value reaches one of the prescribed bounds after op- 
timization. From Table 1, f 6 (1st elastic lead-lag) is an 
active constraint in Case 1. The 3/rev m* and m c (flap- 
ping and torsional moments, respectively) and the thrust 
constraint are also active in Case 1. The remaining con- 
straints stay within or below the prescribed limits. There 
is a significant reduction (11.7 percent) in blade weight, 
from that of the reference blade, in Case 2. Fig. 2 shows 
that the centrifugal stresses in each segment remain below 
the prescribed upper bound, s m ax, and less than reference 
blade stresses for both Cases 1 and 2. 

Table 2 presents a summary of the reference and the 
optimized blade design variables (except for the nonstruc- 
tural masses which are presented in Fig. 3). For Case 1 
all of these design variables except k r remain unchanged, 
from the reference to the optimum. The optimum blade 
therefore remains rectangular in planform. Considerable 
changes, however, occur in Case 2. The increases in EIxx r 
and EIzzr from the reference blade values, in Case 2, are 
caused by the fact that the optimum blade is tapered, with 
X = 1.24, whereas the reference blade has a rectangular 
planform. Therefore, to account for the loss of stiffnesses 
towards the blade tip and to satisfy the frequency con- 
straints, the stiffness values must increase at the root. It 
is interesting to note from Table 2 that the value of the 
design variable GJr remains unchanged after optimization 
although the 3/rev torsional moment is constrained. This 
is because only a rigid body torsion mode is included in 


the analysis. In Case 2 there is a significant reducdon in 
the blade root chord (37.8 percent) resulting in a reduced 
blade solidity (46.7 percent) from the reference blade 
value; the optimum blade in Case 2 with a smaller 
solidity permits less thrust to achieve the same Ct/ct goal 
(Table 1). 

The nonstructural weight distributions of the refer- 
ence and the optimized blades are presented in Fig. 3. 
Reductions are shown in these weights, from reference to 
optimum, primarily over the blade inboard region. The 
weights increase towards the blade tip due to the autorota- 
tional inertia constraint. The inboard reductions of these 
weights are less in Case 1 than in Case 2. Table 2 and 
Fig. 4 present a comparison of the amplitudes of the 4/rev 
vertical and the 3/rev inplane shears, which are the objec- 
tive functions. The 4/rev vertical hub shear is reduced by 
10.9 percent and the 3/rev inplane shear is reduced by 4.1 
percent, from the reference blade values, in Case 1. The 
reductions are naturally larger in Case 2, 89.8 percent and 
55 percent reductions in the vertical and inplane, respec- 
tively. This is due to the fact that this rotor has much 
lower thrust. 

It was of interest to investigate how much of the 
reduction in the vibratory forces and moments were 
caused by reductions in the blade section airloads. There- 
fore, the normal and the inplane components of the total 
aerodynamic forces on the section, resolved with respect 
to the hub plane, and denoted F z and F x respectively are 
plotted in Figs. 5-12. Figure 5 presents the radial dis- 
tribution and Figs. 6 - 8 present the azimuthal distribu- 
tions of F z . These figures show that, for Case 2, op- 
timization significantly reduces the amplitude of F z dis- 
tribution of the advancing blade along the blade radius 
(Fig. 5). The amplitudes of the same distribution around 
the azimuth, plotted at three critical radial stations, 25% 
radius (y/R = 0.25; Fig. 6), the thrust-weighted equivalent 
chord point (y/R = 0.75; Fig. 7) and blade tip (y/R = 0.99; 
Fig. 8), for Case 2, are also reduced. No visible changes 
are observed in the above distributions for Case 1. The 
large reductions in the amplitude of ¥ z distribution in 
Case 2 can be attributed to the lower thrust carried by the 
blade. Therefore, the larger reduction in the 4/rev vertical 
shear in Case 2 is brought about largely by the reduction 
in F z . Fig. 9 presents the radial distribution and Figs. 10 
- 12 present the azimuthal distributions of F x (at the same 
three radial locations as in F z . Only Case 2 shows reduc- 
tions in the amplitudes of F x distribudons. The reason is 
that, in Case 2, the blade is being trimmed to produce the 
same Ct/ct and Cx/g as the reference blade and no con- 
straint is imposed on the thrust and the propulsive force. 
The smaller rotor eventually carries both less thrust and 
less propulsive force. Therefore, F x having a component 
of the thrust as well as the propulsive force embedded in 
it reduces significantly. This study clearly shows that in 
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Case 2 the vibratory root shears and moments are reduced 
partly by proper blade tuning (through nonstructural 
weights) and frequency placements, but largely due to the 
reductions in the blade airloads. The mechanism is dif- 
ferent in Case 1, where blade aerodynamic consideration 
is strongly coupled through the thrust constraint. In this 
case, the mean loads remain unchanged, however, the 
vibratory forces are reduced mainly through proper fre- 
quency placements and blade tuning. 

Effect on Thrust Constraint 

Since the thrust constraint was found to play an im- 
portant role in the optimum blade design, it was of inter- 
est to study its role in greater details. This is done by 
varying the magnitude of the prescribed thrust (by chang- 
ing the value of *n\ Eqn. 4) and comparing the forces and 
moments at the blade root. The results of this study are 
shown in Figs. 13 and 14. Figure 13 presents the varia- 
tions of the root shear forces and Fig. 14 presents the 
variations of the root moments with changes in the thrust 
constraint. It is interesting to note, from these figures, 
that the thrust constraint does not affect the values of the 
forces and moments for values of T ^ 0.54 T re f The mag- 
nitudes of the shear forces and moments increase with the 
increase in the magnitude of the prescribed thrust, as ex- 
pected, above that value. For example there is a near- 
linear increase in the values of the shear forces as indi- 
cated by Fig. 13. It is also of interest to note from Fig. 
13 that the 3/rev inplane and radial shears vary almost 
identically with n. Figure 14 shows gradual increases in 
the values of the 4/rev lagging and the 3/rev torsion mo- 
ment above n = 0.54. The increase in the 3/rev flapping 
moment value tends to flatten out above n = 0.90. In both 
Figs. 13 and 14 the value of n = 1.0 corresponds to the 
100% thrust case (Case 1). 

Effect on Rotor Performance 

Although rotor performance was not considered while 
designing the optimum blades, it was of interest to inves- 
tigate the effect of optimization oh the performance. An 
important performance criterion is the total power re- 
quired for a given task, which is a measure of economic 
efficiency. Therefore, it was of interest to see how it is 
affected. Recalling the fact that the optimum rotor in Case 
2 is significantly smaller than the reference rotor, in Fig. 
15 the power coefficient per unit solidity is compared. 
Fig. 15 indicates that the dynamic tuning did not increase 
the normalized power in Case 1. In Case 2 the solidity 
flexibility reduced the power as the thrust was reduced, 
however, the power loading (per unit solidiy) was more 
efficient for the optimum blade (9.2 percent lower than 
the reference blade) as shown in Fig. 15. This is further 
demonstrated by plotting the normalized section power 


coefficient (normalized with respect to its value at the tip 
of the advancing blade) against Mach number, at some 
typical radial locations, in Figs. 16-20. These figures 
indicate that the general nature of the power distributions 
around the blade azimuth are similar for the reference and 
the optimum blades of both Cases 1 and 2. For example, 
at the thrust- weighted equivalent chord point (Fig. 18) the 
power is evenly distributed between the advancing and 
the retreating blades whereas towards blade tip (y/R = 
0.99; Fig. 20) most of the power is being used by the 
advancing blade. However, there is a definite reduction in 
the power requirement, from reference to optimum, as in- 
dicated by the reduced sizes of the diagrams. The reduc- 
tion being more significant in Case 2. The reduction in 
the power is further investigated through a study of the 
blade angle of attack distributions. Figures 21-24 show 
the angle of attack plotted against Mach number, at the 
same radial locations as the section power coefficient dis- 
tribution. The (L/D)max shown in these Figures, cor- 
responds to the RC - 410 airfoils (Ref. 16) used in the 
reference blade. The figures indicate that the distributions 
are similar between the reference blade and the optimum 
blade in Case 1. However, for Case 2, both figures show 
that the advancing blade moves further away whereas the 
retreating blade moves closer to the (L/D)max. This indi- 
cates that the retreating blade of the optimum rotor is 
operating more efficiently than the reference rotor, caused 
by a shift and increase in the mean section lift coefficient 
on the retreating blade after optimization. This explains 
the efficient power loading in the optimum blade for Case 
2. A similar phenomenon is reported in Ref. 17. 

Effect of Move Limit 

It was observed, during the optimization process 
using linear Taylor series approximation, that the move 
limits (ML) used in the approximate analyses of the ob- 
jective functions and the constraints affected the optimum 
values to a great extent. Depending upon the value of 
ML, several different local minima resulted. One such 
occurrence is shown in Fig. 12, where the objective func- 
tions with a thrust constraint of T> 0.75 T re f (n = 0.75 in 
Eqn. 4) are plotted for different values of ML. The fig- 
ure shows a large increase in the optimum value of the 
4/rev vertical shear, f z , with a decrease in the value of 
ML from 0.1 to 0.02. Similar changes are also in the 
optimum value of the 3/rev inplane shear, f x . However, f z 
is found to be more sensitive to ML than f x . For ex- 
ample, the optimum value of f z reduces by more than 25 
percent as the value of ML is increased from 0.02 to 0.1. 
The corresponding change is around 8 percent for f x . 
This shows that the objective functions used are highly 
nonlinear with respect to the design variables with higher 
nonlinearity in case of the 4/rev vertical shear. Therefore, 
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such nonlinear objective functions and constraints needed 
to be examined more carefully and a study was done on 
the most effective use of ML. In general, there is a trade- 
off between a good local minimum and increased com- 
puter time. For instance, a small move limit provides a 
smoother convergence, but the rate of convergence is 
very slow and often the chances are that it does not con- 
verge to a good local minimum (e.g.. the f z value in Fig. 
25 with ML = 0.02). A large move limit gives a faster 
rate of convergence with more chances of oscillation (in 
case of a nonlinear function), often running into mathe- 
matical singularity (i. e., negative objective function). It 
was therefore decided to use a variable move limit proce- 
dure in which the optimization procedure is started with a 
large ML and the process continued until the approxima- 
tion yields negative objective functions. At this point a 
switch is made to several sets of smaller values of ML. 
The ML value that yields the best local minimum is 
selected to be the optimum design. In the present 
analysis, move limits in the range of 0.1 - 0.01 (10 - 1 
percent) have been used. 

Concluding Remarks 

In this paper an integrated acrodynamic/dynamic op- 
timization approach has been described for rotor blades to 
obtain reduced vibratory root loads. The 4/rev vertical 
shear and the 3/rev inplane shear forces are minimized for 
an advanced articulated rotor blade under forward flight 
conditions. The ’Global Criteria Approach’ is used for 
formulating the multi-objective optimization problem. 
Constraints are imposed on the remaining vibratory loads 
- the 3/rev radial shear, the 3/rev flapping and torsional 
moments and the 4/rev lagging moment. Constraints are 
also imposed on blade natural frequencies, autorotational 
inertia, centrifugal stress and rotor thrust. The program 
CAMRAD is used for the aerodynamic and dynamic 
analysis of the blade and the program CONMIN, along 
with an approximate analysis technique, is used for the 
optimization. The blade is trimmed at each step of design 
optimization. The optimum designs are compared with a 
reference blade. Two optimum design cases were 
studied. In Case 1 the thrust of the optimum rotor was 
constrained to be the same as the reference rotor and in 
Case 2 no constraint was imposed on the thrust. The ef- 
fect of thrust constraint was studied in detail. Rotor per- 
fomance, a criterion, not included while designing the op- 
timum blades, was also studied. The effect of move limit 
used in the approximate analysis on the optimum solution, 
was investigated. 

Following are the conclusions that are made from this 
study: 

1. The Global Criteria Approach is effective in solv- 
ing the multiple objective problem The integrated op- 


timization scheme is very efficient and results are ob- 
tained in 10-15 cycles. 

2. The procedure yields significant reductions in the 
objective functions. The 4/rev vertical shear is reduced 
by 10.9 and 89.8 percent and the 3/rev inplanc shear is 
reduced by 4.1 and 55 percent in Cases 1 and 2, respec- 
tively. The reductions are larger in Case 2 due to the 
lower thrust of the rotor. 

3. The amplitudes of the local section airloads (nor- 
mal and inplane) arc substantially reduced in Case 2, 
proving it to be the driver behind the large reductions in 
the vibratory root forces and moments in this design. 
This airload reduction comes about along with a reduction 
in rotor thrust. In the practical application of a Case 1 
deign, no reduction occurs in the airloads due to the thrust 
constraint and the vibratory forces and moments arc 
reduced through blade tuning and frequency placement. 

4. In Case 1 the optimum blade retains the same nor- 
malized power requirement as that of the reference blade 
and has a more efficient power loading than the reference 
blade for Case 2. 

5. A proper choice of the move limit is important in 
arriving at a better local minimum. 
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Table 1 Contraint Values 




Reference 

Prescribed Bounds 

Optimum 

blade 

lower 

upper 

(Case l) a 

(Case 2) b 

f3 (per rev) 

(flap) 

3.07 

3.05 

3.50 

3.13 

3.31 

f 4 (per rev) 

(flap) 

6.76 

6.50 

6.90 

6.87 

6.90 

f 5 (per rev) 

(flap) 

9.28 

9.25 

9.50 

9.38 

9.25 

f6 (per rev) 

(lead- lag) 

12.63 

12.50 

12.75 

12.75 

12.71 

AI (lb-ft 2 ) 


19.75 

19.75 

- 

20.30 

19.75 

W (lb) 


3.41 

- 

3.41 

3.39 

3.01 

3/rev f r (lb) 


2.81 

- 

2.81 

2.75 

1.08 

3/rev m* (lb-ft) 


0.59 

* 

0.59 

0.59 

0.26 

3/rev m c (lb-ft) 


0.23 

- 

0.23 

0.23 

0.05 

4/rev m z (lb-ft) 


0.70 

- 

0.70 

0.64 

0.18 

Thrust, T (lb) 


298.70 

298. 70 c 


298.70 

158.50 


*Case 1 : optimum design with thrust constraint, T > T rc f 
b Case 2 : optimum design without thrust constraint 
c For Case 1 only 

Table 2 Optimiation results 



Reference 

Optimum 

Optimum 


blade 

(Case l) a 

(Case 2) b 

Elxxr (lb-ft 2 ) 

10277.00 

10277.00 

10300.00 

EIn r ( lb-ft 2 ) 

354.00 

354.00 

378.00 

GJ r (lb-ft ) 

261.00 

261.00 

261.00 

kr (ft) 

0.27 

0.16 

0.17 

X 

1.00 

1.00 

1.24 

Cr (ft) 

0.45 

0.45 

0.28 

4/rev ft (lb) 

0.16 

0.14 

0.016 

Percent reduc- 
tion 0 


10.90 

89.80 

3/rev yi (lb) 

3.17 

3.04 

1.43 

Percent reduc- 
tion 0 

- 

4.10 

55.00 


a Case 1: Optimum design with thrust constraint, 

T 5 Tref 

b Case 2: Optimum design without thrust constraint 
c From reference blade 
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Fig. 1 Rotor blade model. 


CZ3 Reference 



Optimum (Case 1) 

Optimum (Case 2) 1 3.5 


-i 3.0 


4/rev 
vertical shear (f z ) 
(lb) 


- 2.5 



3/rev 
inplane shear (f x ) 
(lb) 


Reference 

Optimum (Case 1) 

x iq 7 Optimum (Case 2) 



Fig. 2 Centrifugal stress distributions. 


Fig. 4 Comparison of objective functions. 


Reference 
Optimum (Case 1) 



Weight, 

w 

Ob) 


Reference 
Optimum (Case 1) 



Fig. 3 Nonsnuctural weight disnibutions. 


Fig. 3 Radial distributions of normal component of 
airload. F t \ ¥ = 90". p. = 0.3 


— Reference 

— Optimum (Case 1) 



Fig. 6 Azimuthal distributions of normal component of 
airload , Ft, = 0.25, (i=0.3. 
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Reference 
Optimum (Cam 1) 



Reference 
Optimum (Case 1) 



Fig. 7 Azimuthal distributions of normal component of 
airload, F z \ y/R = 0.75, ji=0.3. 

— o — Reference 
— a — Optimum (Case 1) 



Fig. 10 Azimuthal distributions of inplane component 
of airload, F x \ VR ~ 0.25, \i=03. 

— o — Reference 
— a — Optimum (Case 1) 
— □ — Optimum (Case 2) 



Fig. 8 Azimuthal distributions of normal component of Fig. 1 1 Azimuthal distributions of inplane component 

airload, F z \ y/R = 0.99, p=0.3. of airload, F x \ y/R = 0.75, }i=0.3. 


Reference 


Optimunv(Case 1) 


Inplane 

component, 

Fx 

(Ib/tt) 



Reference 
Optimum (Case 1) 



Fig. 9 Radial distributions of inplane component of 
airload, F x \ 90*, *1=0.3. 
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Fig. 12 Azimuthal distributions of inplane component 
of airload, F x \ y/R - 0.99, ji=0.3. 



■ — fx (3/rev inplane) 



h 

(lb) 


Reference 
Optimum {Case 1) 



Fig. 13 Effect of thrust constraint on root shear forces. Fig. 16 Variation of normalized local power coefficient 

with Mach number; y/R = 0.25. 


• — m x (3/rev flapping) 



.3r 


.2 - 

Normalized local 
power coefficient, - 1 ' 
c p/ c p,(V=90") 

0 


Reference 

Optimum (Case 1) 

Optimum (Case 2) 



-.1 


0 


.10 .20 .30 
Mach number 


.40 


.50 


Fig. 14 Effect of thrust constraint on root moments. 


Fig, 17 Variation of normalized local power coefficient 
with Mach number; y/R = 0.45. 


.0040 


Normalized 

total power 0035 

coefficient, 


Cp/O 


.0030 


I 

! 


i 

I 


Reference 


Optimum 
(Case 1) 


Optimum 
(Case 2) 


Reference 



Mach number 


Fig. 15 Comparison of normalized total power Fig. 18 Variation of normalized local power coefficient 

coefficient. with Mach number, y/tt = 0.75. 
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Reference 




Fig. 19 Variation of normalized local power coefficient Fig. 22 Variation of local angle of attack with Mach 

with Mach number; y/R- 0.91. number; y/fe = 0.45. 



Fig. 20 Variation of normalized local power coefficient Fig. 23 Variation of local angle of attack with Mach 

with Mach number; y/R = 0.99 number. y/R = 0.75. 




Fig. 21 Variation of local angle of attack with Mach Fig. 24 Variation of local angle of attack with Mach 

number; y/R = 0,25. number; y/R = 0.99 
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EZ3 ML = 0.02 
E2 ML = 0.05 



Fig. 25 Effect of move limit on objective functions. 
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